Microscopic model for transitions from Mott to spin-Peierls insulator in TiOCl 
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On the basis of ab initio density functional theory (DFT) calculations, we derive the underlying 
microscopic model Hamiltonian for TiOCl, a unique system that shows two consecutive phase transi- 
tions from a Mott insulator to a spin-Peierls insulator through a structurally incommensurate phase. 
We show with our model that the presence of magnetic frustration in TiOCl leads to a competition 
with the spin-Peierls distortion, which results in the novel incommensurate phase. In addition, our 
calculations indicate that the spin-Peierls state is triggered by adiabatic phonons, which is essential 
for understanding the nature of the phase transition. 



I. INTRODUCTION 

Mott insulators are materials which owe their insulat- 
ing properties to strong electron correlations 1 . In prac- 
tice, the Mott state is unstable to residual interactions 2 
and in one dimension the dominant instability is the dou- 
bling of the unit cell via the coupling of spins to phonons, 
so called the spin-Peierls transition 3 . 

With the discovery in the early nineties of the first 
inorganic spin-Peierls material, CuGeOj 4 -, followed by 
NaV20g£ a large amount of work£ was devoted to the 
understanding of the complexity of this instability when 
additional degrees of freedom like next nearest neigh- 
bor interactions (CuGeOa) or charge disproportionation 
(NaX^Os) are involved. Moreover, commensurate to 
structurally incommensurate transitions were observed in 
TTF-CuBDTl and CuGeOg&S, when a strong magnetic 
field is applied. 

Recently, the anomalous behavior of TiOC l 10 ' 11 ' 12 - 
the third inorganic spin-Peierls compound with highest 
transition temperature and largest lattice distortion yet - 
has provoked a new debate. Upon cooling, TiOCl under- 
goes two consecutive phase transitions at T C2 = 92 K and 
T Cl = 66 K of second and first order, respectively. These 
two temperatures separate a paramagnetic Mott insula- 
tor at high temperatures (T > T C2 ) from a non-magnetic 
dimerized spin-Peierls state at low temperatures (T < 
T Cl ). In the intermediate region T Cl < T < T C2 , TiOCl 
shows structural incommensurabilities 13 ' 14 ' 15 ' 16 . The 
observation of such a phase without the need of applying 
a magnetic field is unprecedented, making the study of 
TiOCl compelling since this phase must be a result of 
the subtle competition between all interactions present 
in the material. 

TiOCl is also unconventional in the size of the dis- 
tortion rj = 0.18 A observed between neighboring Ti 
atoms along 6 13 ' in the dimerized phase in contrast 
to e.g. polyacetylene where rj — 0.08 A 1 -. Moreover, 
the characteristic frequency of the spin-Peierls phonon 
f2o — 27 meV is comparable to the spin excitation gap 
A ~ 21 meV associated with the static dimerization in 
TiOCl for T < T Cl i&. This indicates that the anoma- 
lous spin-Peierls transition might be at the borderline 



between adiabatic (the phonons behave classically) and 
non-adiabatic (the quantum nature of the phonons is im- 
portant for the phase transition) behavior 18 i 19 i 20 , More- 
over, increasing interest in applying pressure to TiOCl 
21,22,23,24,25 re q U i res decent understanding of the phase 
transitions at ambient pressure. 

Looking at the structure of TiOCl, this system con- 
sists of double layers of titanium (Ti) and oxygen (O) 
in the ab plane separated by chlorine (CI) layers along 
the c direction i&I&Ii (see Fig. Q] (a)). The Ti atoms 
form chains along b with intrachain coupling J\, whose 
projection on the ab plane gives a triangular pattern 
(Fig. [T] (c) and (d)). It has been suggested in the 
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that possible frustrating interchain interac- 



tions {J a within one layer and J c between layers) could 
be responsible for the existence of the intermediate struc- 
turally incommensurate phase with a temperature depen- 
dent modulation wave vector q = ((ft, 0.5 + q%, O) 13 ' 14 ' 16 , 
with < (ft, <72 <S 0.5. But the mechanism driving the 
system to such a phase has not been demonstrated. 

In this paper, we will demonstrate that the behavior 
of TiOCl can be described by a two-dimensional frus- 
trated spin-Peierls model with large magnetoelastic cou- 
pling. We show that magnetic frustration^ is respon- 
sible for the existence of a structurally incommensurate 
phase at T C1 < T < T C2 . Our results indicate that 
the spin-Peierls state is triggered by adiabatic phonons. 
The model parameters in our proposed model are calcu- 
lated within ab initio density functional theory (DFT) by 
the Car-Parrinello projector augmented wave (CP-PAW) 
method^. 

The paper is organized as follows: Section IITI describes 
the model we proposed and the details about how we 
determined the model parameters by DFT methods. In 
Section lLTIl we present our results of the model parameters 
and discuss the consistencies between our results derived 
from the proposed model and the experimental obser- 
vations. The possible mechanism for the structural in- 
commensurability without external magnetic field is also 
discussed and the nature of the phonons is determined. 
Finally in Section HVl we present our conclusions. 
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FIG. 1: (Color online) TiOCl structure and exchange inter- 
action diagrams, (a) View along b direction of TiOCl. Here, 
Ti is shown by blue ball (largest size), O red (smallest size) 
and CI yellow (intermediate size), (b)-(d) Schematic view of 
the structure, showing only Ti ions, (b) Projection on the ac 
plane. Blue (Black)/light blue (light grey) circles denote Ti 
ions in upper/lower layer. The unit cell is shown as a shaded 
region, (c), (d) Projection of the bilayer onto the ab plane. 
Ti chains run along 6. Note the triangular arrangement of 
Ti ions. Couplings between Ti ions are shown as blue (verti- 
cal) (Jb), magenta (diagonal) (J c ) and grey lines (horizontal) 
(Jo)- The various line thicknesses signal the strengths of Ti 
displacements for T < T C1 (d) which are absent for T > T c , 
(c). 



II. MODEL AND METHOD 

The most general two-dimensional frustrated spin- 
Peierls model for TiOCl can be written as follows: 



H — H a + Hf, + H c 



(1) 
(2) 



Hb — Jb ^ (1 + ctbd£,ij,ij+b) Sjj • Si t j + b 
id 

+£^fe + >) 2 +£^ (3) 



H C = J C J2{ 1 + a c d^ ijti+§j± kj S id • S, i+f J± | 

id 



where d£ij t ki denote the displacements of the Ti ions 
in the ab plane d^M = £ij,kl ~ ^%,kv Here ' ^HM = 
— rki \ is the distance between two Ti ions and = 
r° + Uij is the Ti position for T < T Cl (Fig. [T] (d)) ob- 
tained by adding a displacement Uy to the original Ti 
position for T > T C2 (Fig. [T] (c)). ab and a c denote 
the spin-phonon couplings and K is the elastic constant. 
Mb and M c are the effective masses obtained by weighted 
summation of the masses of Ti, O, and CI atoms and Pb,ij 
and P c ,ij denote the momentum at site along b and 
c respectively. Displacements along a are not included 
since they are absent in the low-symmetry structur o 13 ' . 

In order to obtain the model parameters from first 
principles analysis, a fully relaxed T > T C2 crystal 
structure is needed. Calculations with various func- 
tionals (generalized gradient approximation (GGA) and 
GGA+U) with the CP-PAW method showed that the 
lattice parameters as well as distances and angles from 
GGA+U are more consistent with experimental results 
than those from GGA, demonstrating the importance of 
including electron correlation effects for the description 
of the lattice properties of TiOC1^2^, even if it is done at 
the simple level of GGA+U. Moreover, our CP-PAW cal- 
culations are all carefully checked to be fully converged 
with respect to the size of the basis set. 

The magnetic exchange coupling constants Ji = 
Jb,J a , Jc were obtained from total CP-PAW energy dif- 
ferences between the relevant ferromagnetic (FM) and 
antiferromagnetic (AFM) Ti spin configurations -Efm — 
-E^afm by relating them to the Ji through mapping to 
the classical Heisenberg model. The determination of 
the spin-phonon couplings ab and a c was more elab- 
orate. Out of the phonon spectrum of TiOCl for the 
high-symmetry (T > T C2 ) and the low-symmetry (dimer- 
ized) (T < T Cl ) phased, the spin-Peierls phonon was 
identified and found to be degenerate in energy (A fl 
and B„). The displacement vector for the Ti atoms 
u , = (0.0, +0.0223, ±0.0403)<5, corresponding to the 
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FIG. 2: (Color online) Exact diagonalization results for 
spin susceptibility, (a) Susceptibility, calculated with the 
model parameters determined at U = 1.65 eV. The 1-, 2- 
, and 3-chain geometries are shown in the inset. Verti- 
cal/diagonal/horizontal bonds denote the Jb/Jc/Ja exchange 
coupling, respectively. Blue (Black)/light blue (light grey) 
circles denote Ti ions in upper/lower layer as Fig. [TJ (b) Sus- 
ceptibility calculated with the model parameters of Ref. [32| 
which clearly disagrees with the susceptibility calculated for 
the experimentally determined coupling Jt, = 660 K. 



spin-Peierls phonon, was then calculated by diagonaliz- 
ing the dynamical matrix, where 8 is the order parameter, 
with 8 = 2 corresponding to the displacements of the Ti 
atoms in the dimerized phase. We obtained a c by cal- 
culating the difference in energy between the FM and a 
spin S — 2 AFM spin arrangement of the unit cell dou- 
bled in both b and c directions and distorted according to 
the displacement vectors Uy- for 8 = 2. ab was calculated 
analogously. We performed the calculations for various 
U values and found that J c and J a are always FM and 
relatively small compared to the exchange coupling Jf,; 
the latter is always AFM and decreases monotonously 
as U increases, roughly fulfilling the relation J ~ 1/U. 
We have also performed LAPW calculations with the 
WIEN2k 30 code to check these results and we obtain 
good agreement between both methods. 

The elastic constant K was calculated by the expres- 
sion E{S)=E (0) + {).{){)72NK8 2 /2 where N is the num- 
ber of Ti ions in the unit cell and the numerical factor in 
the 8 2 term arises from the relationship between 8 and 
d£ (see equation UJ. E (8) is the CP-PAW ground state 
energy obtained for distorted lattices according to the 
spin-Peierls vectors given above. K was found to be 
almost independent of the choice of U. 



III. RESULTS AND DISCUSSIONS 

For the moderate value U = 1.65 eV— , we obtain 
the model parameters Jb = 660.1 K, J c = —16.7 K, 
J a = -10.5 K, a b = 1.73 A" 1 , a c = 8.86 A" 1 and 
K = 0.292 eV/A 2 . The reliability of these parameters 
is demonstrated by comparing the spin susceptibilities 
among single-, double- and triple-chain Heisenberg mod- 
els calculated by exact diagonalization (see Fig. [3] (a)). 
The behavior of the spin susceptibility at moderate to 
high temperatures remains almost unchanged when inter- 
chain couplings J c and J a are gradually introduced from 
single- to triple-chain models. This insensitivity of the 
susceptibility to the interchain couplings explains why 
initially TiOCl was described as a good realization of a 
one-dimensional spin 1 /2 system^ since the experimen- 
tal data could be perfectly fitted to a one-dimensional 
spin 1/2 Heisenberg chain. Our results show that the 
experimental susceptibility is well described by a frus- 
trated model. In contrast, a recent estimate of the mag- 
netic interactions given by Macovez et al. 32 for TiOCl 
with J b = 353 K, J a = -7.4 K and J c = 40.8 K shows 
in our ED calculations considerable deviations from the 
experimental susceptibility as shown in Fig. [5] (b). 

The obtained values of J a and J c are the key for under- 
standing the two consecutive phase transitions in TiOCl. 
The cooperation of J a and J c (2J C + J a = — 43.9K) in 
competition with Jb on a Ti site is in the range of ener- 
gies where the incommensurability happens in the system 
(between T Cl = 66 K and T Cl = 92 K). Starting from the 
non-magnetic dimerized phase, this ferromagnetic inter- 
action together with thermal fluctuations suppresses the 
formation of dimer states in favor of formation of spin 1/2 
solitons leading to a structurally incommensurate phase. 
In Fig. [3] we present an illustration of the process. As 
temperature increases, the ferromagnetic coupling breaks 
the central spin zero dimer into a triplet state, and the 
neighboring spins get partially polarized. These spins 
feel the attraction (repulsion) of the central Ti S = 1/2 
spins depending on whether their interaction is favorable 
(unfavorable) according to J a , J c and Jb- Consequently, 
helped by the thermal fluctuations, the neighboring Ti 
atoms will displace towards (away from) the central spins 
(see Fig. [3](b)). This displacement will propagate to the 
next nearest neighbor atoms so that a sinusoidal atomic 
displacement (soliton) over many unit cells along the a 
and the b direction will appear, which defines the incom- 
mensurability along a and b observed experimentally. 

The couplings in TiOCl reveal still further important 
features of the spin-Peierls phase transition. Since Jb 3> 
J a , J c , we can define a dimensionless magnetoelastic cou- 
pling constant, which uniquely determines the properties 
of the system, as A = Jba 2 /K with Adft = 0.58 as calcu- 
lated from our model parameters for TiOCl. This value 
is significantly larger than A = 0.32 for polyacetyleneii. 
With A we can now tackle the question whether the spin- 
Peierls transition is adiabatic or non-adiabatic. 

We analyze Hb (equation ([3]) ) in the adiabatic limit by 
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FIG. 3: (Color online) Explanation of incommensurate struc- 
ture. Here, blue (black) /light blue (light grey) circles denote 
Ti ions in upper/lower layer as Fig.fj] (a) Dimerized structure 
at temperatures below T Cl = 66 K. The antiferromagnetic in- 
teractions and spin-phonon coupling along b lead to the for- 
mation of stable singlets, (b) Above T C1 = 66 K, thermal 
fluctuations lead to a competition between antiferromagnetic 
interactions along b and frustrating ferromagnetic interactions 
along a and c. We illustrate the effect on the atomic positions 
by showing the result of a spin flip in the central Ti dimer 
(black solid spins). The resulting attractive (inward arrows) 
or repulsive (outward arrows) interactions (gray (light grey) 
line J a and magenta (dark grey) J c ) lead to modulations of 
the Ti atom positions along a and b directions, explaining the 
formation of an incommensurate structure. 



(i) using the mean field (MF) approximation after apply- 
ing a Jordan- Wigner transformation and by (ii) perform- 
ing exact diagonalization (ED) calculations for chains 
with up to 24 spins at T = 0. In Fig. 2] (a), we show the 
distortion of the lattice rj as a function of the coupling 
constant A for the two calculations. The MF approxima- 
tion predicts a Amf ~ 0.56 to induce the experimental 
distortion. Note that this Amf is very close to our ab 
initio determined value Adft- ED results qualitatively 
agree with the MF results and predict a Aed ~ 0.64 
to describe the experimental distortion. This agreement 
between Adft and Amf or Aed demonstrates the high 
degree of consistency in our calculations. In Fig. 0] (b), 
we show the temperature dependence of the phase tran- 
sition within the MF approach for different A values. We 
observe that the spin-Peierls phase transition is of second 
order and the critical temperature T$p lies between 100 
and 230 K depending on the choice of A, which is not far 
away from the experimental T C2 . As is well known, in- 
clusion of interchain couplings and strong correlation in 
the calculation further reduces Tgp to the experimental 
value T C2 . 

We now consider Hf, (equation l[3"])) in the non- 
adiabatic limit. After quantizing the lattice degrees of 
freedom by introducing the phonon creation and annihi- 
lation operators and integrating out the phonon degrees 
of freedom, an effective J\-Ji model is obtained, where 
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FIG. 4: (Color online) Predictive power of calculated model 
parameters, (a) Lattice distortion r\ as a function of the 
magnetoelastic coupling A obtained from various approaches: 
Heisenberg model in mean field approximation (black curve) 
and exact diagonalization calculations for systems up to 24 
sites (symbols) for A = 0.56 (triangles), 0.64 (circles), and 
0.72 (squares). Note that sizes of N = 24 are large enough to 
reach converged results. The horizontal dashed line denotes 
the experimental rj. (b) Lattice distortion r\ as a function of 
temperature for various A obtained for Hi, in the mean field 
approximation. The second order nature of the T C2 phase 
transition is evident from the curves. 



J 2 = (Jba b ) /4 Mfcfi| 0.036 J fc and J x = J b + 2 J 2 w 
1.07Jh (Qq — 2 yj K/Mb is the bare phonon frequency). 
Since the ratio J%l J\ < 0.24-, non-adiabaticity will not 
affect the dimerization in TiOCl (in contrast to the situa- 
tion in CuGe0 ^ 19 i 2 ° ). In fact, with our calculated model 
parameters, we find that Qq/Tsp ~ 2 fulfills the condi- 
tion Qo/Tgp < 2.32£ for the occurrence of a soft phonon. 
Out of these two limiting calculations, we conclude that 
the main features of TiOCl can be well described in the 
adiabatic limit, in agreement with recent experiments^. 



IV. CONCLUSIONS 

The presented work shows that TiOCl is described by 
a frustrated spin-Peierls model on an anisotropic trian- 
gular lattice with intrachain AFM and interchain FM in- 
teractions and large magnetoelastic coupling. This is in 
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contrast with previous suggestions on the magnetic inter- 
actions. With this model we can account for the Mott to 
spin-Peierls insulator phase transition in TiOCl and for 
the existence of intermediate zero-field structural incom- 
mensurabilities. This work also shows that the combina- 
tion of ab initio DFT calculations with effective models is 
a powerful tool to describe the microscopies of this Mott 
insulator. 
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APPENDIX A: MEAN-FIELD MODEL 
CALCULATIONS 

By means of Jordan- Wigner transformation to Hi, 
(equation a spin-less fermion model is obtained, 



Hi 
Jb 



£(l + (-1)* A 



- (c\c l+1 + h.c. 



Uj 

1 2 



n i+ i - 



4-2 A ' 



(Al) 



where A — ctb£,ij.ij+bi cj (c;) is the creation (annihila- 
tion) operator of an electron and rii — c\ci. Applying 
the Hartree-Fock approximation to the interaction part 
and with the help of Fourier transformation and unitary 
transformation, three coupled self-consistent equations 
can be obtained from M- 1 = 0, W- 1 = and ■§£ I = 
where F is the free energy, and the chemical potential /x is 
adjusted to yield the correct filling, i.e., N e = -^^lnS. 



Here, (5 — 1/ksT and 3 is the grand canonical partition 
function 



hi 



S = 2 In [l + e -' 3 (4 F W-p) 



(A2) 



k.± 



where the band structure e^ F (k) = ±\/E (k) and 
£(/.-) = (] - Ar[±-b + 5b 



i I i A) 2 ( i - b - Sb 



(A3) 



1-A 




--b ) - (5b) 2 ) 2 ens/,-. 



bo and 5b is defined as f^c|cj-|_i\ + (c\ +l Ci + 2^j jl and 

r/cjcj+i^ - ^c] +1 Ci +2 \) /2, respectively. Finally, the 
self-consistent equations read 



A = - 



^E{/( £ HF (*))-/(4p(fc))} 

e fc 

1 dE (k) 



ib 5b 



b 



2^/WJk) dA 

J:E{/( e HF(fc))-/(4F (fc))} 

1 dE (fc) 



(A4) 



-4A5b 



N, 



db 

J-E{/( e H F W)-/(4 F (fc))} 

1 dE (k) 
where / (e) is the Fermi distribution function. 
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